Parity projection of QCD sum rules for the nucleon 
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The nucleon and its negative-parity excited states are examined in a maximum entropy method 
analysis of QCD sum rules. First, we rederive the parity projected sum rules for baryons using 
forward correlation functions. Doing this, the method is generalized so that higher order operator 
product expansion (OPE) terms can be calculated unambiguously. We then apply this approach to 
the nucleon channel taking into account all known first order ot s corrections to the Wilson coefficients 
of the OPE. As these corrections have turned out to be large, we suppress them by using a phase- 
rotated Gaussian kernel. Simultaneously, this phase rotation strongly suppresses the continuum 
contribution and improves the convergence of the OPE. The resulting sum rule has the interesting 
feature that it is dominated by the term containing the chiral condensate of dimension 3. Analyzing 
this sum rule by the maximum entropy method, we are able to extract information of both the 
positive and negative parity states. 
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I. INTRODUCTION 

The properties of baryons have been investigated by 
QCD sum rule studies already for more than three 
decades. Historically, this method was developed by Shif- 
man et al. [I], Q and applied to baryonic channels by Ioffe 
0|. Since then, the nucleon sum rules were continuously 
improved by includinghigher orders in the perturbativc 
Wilson coefficients 0t3 or non-perturbative power cor- 
rections 

UBS El- Furthermore, QCD sum rules also 
have been used to investigate the nucleon properties in 
nuclear matter [P2l4l4j or at finite temperature [l5[ . How- 
ever, even in the vacuum, it has not been an easy task to 
appropriately extract information on the nuclconic states 
because of several still not completely resolved issues. 

Specifically, in analyzing the nucleon sum rules, we face 
four major problems. The first problem is that the nu- 
cleon operator couples to both positive and negative par- 
ity states. The contributions of these states can poten- 
tially disturb the analysis and especially lower the relia- 
bility of the extraction of the negative parity states. To 
overcome this difficulty, Jido et al. [HI and Kondo et al. 
[l7l | have constructed certain parity projected sum rules 
which separate the contributions of the positive and neg- 
ative parity states. However, in these analyses, the first 
order a s corrections and the contributions of the non- 
perturbative power corrections above dimension 6 were 
not considered. Furthermore, some technical problems of 
the approach proposed in flq ] were pointed out in [l7| . 
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which are related to a not fully justified use of the for- 
ward correlator instead of the time-ordered one and en- 
suing ambiguities in the treatment of the poles at zero 
energy appearing in the higher order terms of the opera- 
tor product expansion (OPE). 

The second problem is that the a s corrections tend to 
be large [l8[, which seriously put the convergence of the 
Wilson coefficients as an expansion in a s into question. 
This issue is especially serious for the parity projected 
sum rules, for which one cannot avoid the perturbativc 
dimension term, where the a s corrections are large for 
all available local interpolating fields. 

The third problem is the large contribution of the con- 
tinuum, which makes it difficult to investigate the prop- 
erties of the nucleon and its excited states. Indeed, in 
a previous study [Tc| . it was shown that the Borel win- 
dow of the chiral even sum rule is only very narrow if it 
opens at all. Although this difficulty can be avoided by 
only using the chiral odd terms of the OPE, the chiral 
even terms have to be used to separate the positive and 
negative parity contributions. 

The fourth problem is mainly concerned with the ex- 
traction of the lowest negative parity state. Namely, it is 
not known which interpolating field predominantly cou- 
ples to such a resonance and which one merely to con- 
tinuum states. In the traditional analysis, it is necessary 
to model the spectral function by a "polc+ continuum" 
ansatz assuming the existence of a sharp peak, which can 
contradict the actual situation and therefore lead to mis- 
leading results. As we have shown in a previous study 
[l9l ] , the spectral functions corresponding to different in- 
terpolating field operators can show a distinct behavior. 
Some of them do not contain any pole structures. 

To overcome these above mentioned difficulties and to 
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obtain more reliable results, we construct the parity pro- 
jected sum rule including first order a s corrections di- 
rectly from the forward correlation function to resolve 
the problems pointed out in [l7j . and analyze the nu- 
cleon channel by using the maximum entropy method 
(MEM). The MEM analysis of QCD sum rules is a novel 
approach which does not require the "pole+ continuum" 
ansatz and has been applied to the p meson [20| and the 
nucleon channel fl9| in the vacuum and to the charmo- 
nium [2l| and bottomonium (22j channels at zero and 
finite temperature. However, when using conventional 
kernels (Borel or Gaussian), because of the large a s cor- 
rections, the perturbative expansion is not under control 
and it is quite possible that higher order corrections spoil 
the results. Hence, we construct a new sum rule by using 
a phase-rotated Gaussian kernel, in which the a s cor- 
rections are suppressed. Additionally, this kernel also 
suppresses the large contiuum contribution and thus im- 
proves the extraction of the nucleon properties through 
the spectral function. Moreover, this sum rule has the in- 
teresting feature that it is dominated by the chiral con- 
densate term and thus is suitable for investigating the 
relation between the splitting of the positive and nega- 
tive parity states and chiral symmetry. 

The paper is organized as follows. In Sec. [Hi we intro- 
duce the parity projected sum rules, discuss the problem 
of using the forward correlator for constructing the sum 
rules and explain how this issue can be resolved. Then, 
we apply our method to the nucleon sum rule taking into 
account all known first order a s corrections. The max- 
imum entropy method is briefly introduced in Sec. IIII1 
after which in Sec. IIVI the results of the analyses are 
outlined. Summary and conclusions are given in Sec. [V] 



II. PARITY PROJECTED BARYONIC SUM 
RULES 

A. Parity projection of the correlation function 

In QCD sum rules, one usually studies the properties 
of the time ordered correlation function: 
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Here, r\ is an interpolating field constructed from quark 
and gluon operators to have the same quantum numbers 
as the hadron of interest. Making use of the analyticity 
of this correlation function, we can obtain the dispersion 
relation given below: 
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where p(t) denotes the hadronic spectral function, which 
contains contributions from all possible physical states 
coupling to r). Note that we ignore here possible subtrac- 
tion terms, which play no role in the following discussion. 



Eq.([2]) is the usual starting point for QCD sum rule anal- 
yses. 

In baryonic channels, the interpolating field rj carry- 
ing positive parity couples to both positive and negative 
parity baryon states, \n + (q)} and \n~(q)): 



(0\ V (x)\n+( q )) 
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Here, u±(q) are Dirac spinors of positive and negative 
parity and correspond to the strength of the coupling 
of j] to the state \n). Hence, the correlator of 7/ will have 
the following form: 



U(q)=i / d*xe iqx (0\T[r}(x)rj(0)]\0) 
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Here, p + (_)(m) contains the contributions of only pos- 
itive (negative) parity states. When only the sum rule 
constructed from either Ili(q 2 ) or n 2 (g 2 ) is used, it is not 
possible to determine the parity of any structure appear- 
ing in the respective spectral function without additional 
assumptions, as both positive and negative parity states 
contribute to it. Thus, the problem of parity projection 
boils down to consistently disentangling the contributions 
of positive and negative parity to Hi(q 2 ) and n 2 (g 2 ). 

To separate each parity contribution, the forward cor- 
relator in the rest frame (q = 0) was defined in [l6j . 



n for (g )=i J d 4 xe lllx 9(xo)(0\T[ V (x)rj(0)}\0} 

= l0 vl OT M + n 2 OT (<zo), 



9=0 (5) 



where the essential difference to Eq.Q is the insertion 
of the Hcaviside step- function 9(xq) before carrying out 
the Fourier transform. Note that in [l6| . this expression 
was called the "old fashioned" correlator. In the present 
paper, we prefer to call it "forward correlator". 

Using Eq^.Q and the Heaviside step- function, it can 
be shown that the functions nf^go) and n 2 or (go) con- 
tain only poles in the positive go region and are analytic 
for Imgo > 0. The details of how this can be done are 
explained in Appendix A. Furthermore, by applying the 
projection operator ^(70 ± 1) = P ± to the forward cor- 
relator of Eq.© and taking the trace over the spinor 
indices, we are able to construct functions that contain 
only positive or negative parity states, as 

lTr[P±n for (go)] = n^ or (go ) ± H 2 OT (g ) = n±(g ) 

! (6) 

p±(m) dm. 

g - m + ie 
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The imaginary parts of W* 1 (qo ) defined above then give 
the desired parity projected spectral functions: 



-ImFfgo) 
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p±(qo)- 



(7) 



These can, however, not be calculated directly because 
perturbation theory is not reliable in the low qo region, 
but only at |qo| — > oo. 



B. Construction of the sum rules 

Next, we will discuss the construction of the parity 
projected sum rules. For this purpose, we define two ex- 
pressions describing the correlation function (qo)- The 
first one is obtained from the OPE at high energy and is 
thus called nQ PE (qo) in what follows. The second one, 
denoted as IIp hys (qo), is expressed using the hadronic 
degrees of freedom in the physical energy region of qo- 

To construct the sum rules, we consider the contour 
integral 



dqo[ll 
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Xqo)]W(q Q ) = 0, (8) 



where the contour C is given in Fig. [TJ W(qo) must be 
analytic in the upper half of the imaginary plane and real 
on the real axis. As long as it satisfies these conditions, 
it can be chosen arbitrarily. That Eq. (J8J) gives follows 
from the analyticity of both nQ PE (q ) and IIp hys (qo). 
In other words, there are no poles or cuts inside of the 
contour C. Now, we know from asymptotic freedom 
that the perturbative expression nQ PE (q ) approaches 
the hadronic one at |qo| — > oo. Thus, the integrand of 
Eq.© vanishes on the half circle of C, whose radius is 
taken to infinity. Therefore, we are left with the section 
of C which runs along the real axis and can write down 
the sum rule as 



dqoTlo PE (qo)W(q ) 



dqo^ hys .(qo)W(q ). 



(?) 

Finally, taking the imaginary part of the above equation 
we get 



dq Q ^Im[nop E (qo)]W(qo) 

dqoPp hys .(qo)w(qo), 



(10) 



where we have made use of the definition Pp hys (qo) = 
i±mrip hys (q ). Moreover, on the right-hand side, we 
have exploited the fact that IIp hys (qo) only has poles 
on the positive side of the real axis (see Eq.©) and have 
hence restricted the integral to this region. 

Here, we comment on some important differences be- 
tween the sum rules constructed in this paper and the 
one considered in the previous study of |16| . In contrast 




FIG. 1: The contour integral C on the complex plane of the 
variable qo, used in Eq.©. For the actual calculations, the 
radius of the half circle of C is taken to infinity. The wavy 
line denotes the non-analytic cut (or poles) of II ± (qo) on the 
positive side of the real axis. Compared to the discussion in 
the main text, we have here slightly shifted the contour away 
from the real axis into the upper half of the imaginary plane 
for better visuality. 



to Eq. tjTO)) . the authors of [16| have also restricted the 
integral on its right-hand side to positive values. Even 
though this procedure is roughly correct, it leads to am- 
biguities in the higher order OPE terms, which in the 
chiral limit have poles at qo = 0. Furthermore, in [f|| 
the imaginary part of the time ordered correlator (with- 
out 9(xo)) was used instead of ImIlQp E (qo), which in 
principle should be derived from the forward correlator 
of Eq.©. As was pointed out in [f^l, it is not entirely 
clear whether this prescription is justifiable. Therefore, 
in this study, we implement two essential improvements 
as compared to : I) We derive all results directly from 
the forward correlator of Eq.©. 2) We do not restrict 
the region of integration of the OPE side of Eq. (fTUf to 
positive values and therefore remove the ambiguities that 
might occur for higher order OPE terms. 

As a last point, we briefly discuss another parity pro- 
jected sum rule which was proposed in 17 1. There, the 



sum rule is constructed from the retarded correlation 
function, whose spectral function has poles corresponding 
to positive (negative ) parity states in the positive (nega- 
tive ) qo region, respectively. Hence, the contributions of 
the positive and negative parity states are separated into 
the different energy regions of the integral in Eq.©. One 
then can disentangle positive and negative parity states 
by using suitable combinations of several different ker- 
nels. In principle, it is also possible to perform an MEM 
analysis for such combined sum rules and we have indeed 
carried out a trial study to do this. However, in the ac- 
tual analysis, it has turned out to be difficult to extract 
information of positive (negative) parity states which is 
not affected by negative (positive) parity states. The 
reason is that the MEM results have the tendency to be 
symmetric with respect to the positive and negative en- 
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ergy contributions, which can cause artificial peaks [231 ] . 
Therefore, to avoid this problem we use the method of 
fl6| in this paper. 



C. The parity projected nucleon sum rule 

In this subsection, we construct the parity projected 
sum rule of the nucleon including all known a s correc- 
tions. For the nucleon, there are two independent local 
interpolating operators: 



2 (u Ta (x)C-/ 5 d b (x))u c (x), 



m (x) = e abc (u Ta (x)Cd"(x)) l5 u c (x). 



(11) 



(12) 



Here, abc are color indices, C stands for the charge con- 
jugation matrix, while the spinor indices are omitted for 
simplicity. A general interpolating field can thus be ex- 
pressed as 



r){x) = ar]i(x) + /3r] 2 (x), 



(13) 



where a and ft are real parameters. From the time or- 
dered correlator 0, [H| , we can obtain the imaginary part 
of the forward correlator of Eq. (J5)) up to dimension 8 in- 
cluding first order a s corrections. The details of this 
calculation are explained in Appendix [A] 

To construct the final sum rule, we have to specify 
the kernel W(qo). For finding its most useful form, we 
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have tested several kinds of W(qo), such as exp(— 
and j=== cxp(— < ~ 9 " 4t ^ ), which correspond to the Borel 

respectively Doing this, we 
the a, correc- 



/4ttt 4t 

and Gaussian sum rules 
however found that for these sum rules 
tions for the chiral even terms are very large and thus 
far from being under control. For example, the relative 
magnitude of the first order correction of the perturba- 
tivc term docs not depend on a and j3 and reaches values 
as large as 90 % of the leading order term, when taking 
a s = 0.5. Furthermore, the Wilson coefficient of the four 
quark condensate (qqqq) of dimension 6 contains a s cor- 
rections that can be even larger than the leading order 
term and seem to be converging only in a narrow re- 
gion around /3/a = —1. These large corrections seriously 
put the validity of the above sum rules into question. 
Moreover, there is still another ambiguity present in the 
dimension 6 term, which comes from the approximation 
of reducing (qqqq) to (qq) 2 K, k being a dimensionless pa- 
rameter. Although this value is known to be 1 in the 
large N c limit p], y] , its value is only weakly constrained 
at A c = 3 [II. 

As we will explain below, it is possible to improve 
the above situation by choosing an appropriate kernel. 
To do this, we follow the method proposed by Ioffe and 
Zyablyuk [2J|, who have constructed a new class of sum 
rules by using the phase rotated complex variable q 2 e 
instead of the real q 2 . One advantage of this approach 



lies in the possibility of suppressing certain terms of the 
OPE by choosing some specific value of 9. 

To apply this method to the parity projected sum rules, 
we change 
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to 



W(q ,9)dq Q 
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dq 
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for which we note that the additional phase factor ap- 
pearing in front of dqo is needed due to the change of 
the integration variable in Eq. (|10[) . Furthermore, let us 
make a comment on the analyticity of this kernel, which 
is essential for the formulation of the sum rules and thus 
might be a matter of concern as only the real part is re- 
tained in Eq. (fT5"|) . However, as go is real for the integrals 
in Eq. (|10[) . one can rewrite the kernel as 



W(q Q ,6) 



1 1 



q e 2ie exjp( 



-2i0 _ g \2 



4r 



q e 2ie exp 
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(16) 



which of course is an analytic function and can thus be 
continued into the imaginary q plane. As one further 
point, note that the kernel of Eq. (fT5"|) can only be used in 
a limited range of 9, as for too large values of the phase, 
contributions from high-energy states to the sum rule 
will no more be exponentially suppressed. Specifically, 
one can easily convince oneself that 9 should be chosen 
as|0|<§._ 

Introducing the phase rotated kernel explained above 
has two essential advantages: 1) By choosing an appro- 
priate value for the phase 9, one can largely reduce the a s 
correction of the perturbative term and its relative con- 
tribution as a whole to the sum rule. 2) Taking W(qo, 9) 
to be an odd function of qo, both the leading and a large 
part of the subleading a s corrections of the dimension 
6 term in the OPE vanishes. The contribution of this 
term to the sum rules is thus strongly suppressed, which 
means that the uncertainty originating from the value of 
the four-quark condensate has only a small influence on 
the results extracted from the sum rules. 

Substituting the Eq . ([A"22|) and Eq.(rj5]) into Eq. ([TUl) . 
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Here, Gf r (s,r) and G5, 01 '(s,t) arc obtained as follows: 
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where, G , G 0a , G 3 ,G 3cts , C 4 , G 5 , G 6 , C 6as and G 7 are 
given as 
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The values of the vacuum condensates appearing in the 
above equations are shown in Table. |U For all the above 
results, the renormalization constant [i has been set to 1 
GeV. Related to this choice, we use the value of a K = 0.5 
which corresponds to the energy scale of 1 GeV [26|]. To 
obtain the concrete form of G ° t (s,t) and G2° r (s,r), we 
have made use the fact that the contributions of the path 



(qq) 




(qga ■ Gq) / (qq) 


— (0.24 ± 0.01) 3 GeV 3 


0.012 ± 0.0036 GeV 4 


0.8 ±0.2 GeV 2 



TABLE I: Values of the parameters appearing in the OPE, 
taken from [25i. 
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FIG. 2: The paths of Pi and P%, which are used for calculating 
the OPE result given in Eqs. (|18l20|) . The wavy line denotes 
the non-analytic cut (or poles) of 11* (go) on the positive side 
of the real axis. 
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FIG. 3: The ratios of a s corrections to leading order terms 



of dimension and 3: C< £"° (9) and (9) at a = 1.0 and 
/3 = —0.9. These values are shown as a function of s, with 
r fixed to 0.5 GeV 4 . The solid and dotted lines shows 

Co 

at 9 — 0.1087T and 9 — 0, respectively. The dashed line gives 

c 0o 



c 3o 



Co 



which does not depend on 9. 



Pi are as same as those of P2 in Fig. [2] due to the Cauchy 
theorem. 

To determine a suitable value of for which the a s 
corrections of the pcrturbative term are reduced, we find 
that || and 20tan(50) in Go Qs cancel each other out at 
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FIG. 4: Gqpe ( s ' r ) an d the respective OPE contributions 
from dimension to 7 at r = 0.5[GeV 4 ], a = 1, y9 = —0.9 and 
6* = 0. The thick solid lines denote Gope( s i t ) an d the thin 
solid (dashed) lines the chiral even (odd) OPE terms. 



around 9 = 0.1087T, while I\ u gives only a small contribu- 
tion. Hence we set the value to 0.1087T throughout this 
work. To illustrate the effect of the introduced phase, 
the ratios of a s corrections to leading terms at dimen- 
sion and 3: %^(6> = 0.108tt), %^(6> = 0) and %^ 
(which does not depend on 9) are shown in Fig. [3] for 
t = 0.5 GcV 4 . As can be observed from this figure, the 
contribution of the a s term of dimension is reduced to 
5 % of the leading term, which shows that the conver- 
gence of the perturbative expansion is significantly im- 
proved. 

The contributions of each term to G ± (s,r) at 6 = 
and 9 = 0.1087T are given in Figs. [4]and[5j respectively. 
It firstly should be noted here that for both cases, the 
contribution of dimension 6 is small and that the uncer- 
tainty of the four-quark condensate therefore does not 
lead to much ambiguities for these sum rules. This is 
caused by our use of and an go-odd kernel, which elimi- 
nated the leading order contribution of dimension 6. For 
9 = 0, the dimension and 3 terms are dominant, which 
means that not only low lying nucleon states but also the 
continuum (which mainly originates from the perturba- 
tive dimension term) largely contributes to G ± (s, r). 

This situation makes it difficult to extract information 
on nucleonic properties. Especially, for the negative par- 
ity states, the extraction is highly complicated due to 
the strong cancellation between the terms of dimension 
and 3. On the other hand, from Fig. [5J it can be 
seen that in the phase rotated sum rule the chiral con- 
densate term clearly gives the dominant contribution to 
the OPE. Therefore, unlike the sum rule at 9 = 0, the 
disturbance due to the contribution of the continuum is 
considrably reduced. In Fig. [51 it can also be observed 
that the difference between the OPE for positive and neg- 
ative parity is mainly caused by the switched sign of the 



FIG. 5: Same as for Fig. H but for 9 = 0.108tt. 

quark condensate term. Since the spontaneous breaking 
of the chiral symmetry is considered to cause the mass 
splitting of the positive and negative parity states of the 
nucleon, this sum rule is likely to be useful for investi- 
gating the relation between the chiral symmetry and the 
mass splitting between positive and negative parity nu- 
clcons, which we plan to investigate in detail in a future 
work. 



III. THE MAXIMUM ENTROPY METHOD 

In this section, we briefly introduce the maximum en- 
tropy method (MEM), and explain how this approach 
applies to the analysis of QCD sum rules. The equation 
to be analyzed can be written down as 

p DO 

Gope(s)= / W{x,qo)p{q )dq . (22) 
Jo 

Here, x stands for the parameters occurring in the kernel 
W(x,go)j such as s and r. In a QCD sum rule analysis, 
one aims to extract information on the spectral function 
from Gope(^)- However, this function is only known 
with limited precision due to the uncertainties of the val- 
ues of the various vacuum condensates. Furthermore, 
only the limited region of x, where the OPE converges, 
can be used for the analysis. Therefore, rigorously solving 
Eg. (|2"2"]l is an ill-posed problem. Nevertheless, the MEM 
technique enables us at least to statistically determine 
the most probable form of p(qo). 

For this purpose, we define P[p\GH], which is the con- 
ditional probability of p given Gope and H, H represent- 
ing additional knowledge on the spectral function such as 
positivity and asymptotic values. Using Bayes' theorem, 
P[p\GH] is rewritten as 

P[G\pH]P[p\H) 
mGH] = P[G\H] ' (23) 
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where P[p\H] is the so-called prior probability, and 
P\G\pH\ stands for the likelihood function. The most 
probable form of p(lo) is obtained, by maximizing 
P[p\GH] as given in Eq. (f2"3]) . For carrying out this task, 
one needs the concrete functional forms of P[p|H] and 
P[G\pH], which can be expressed as [13, H3 



P[G\pH] cx e~ L[p] , 
P[p\H] oc e aS[pl . 



(24) 



Here, a is a real positive number, which determines the 
relative weight between L[p] and S[p\. The functional 
L[p] is normally used for \ 2 fitting and can be given as 



L[p] 



1 



dx x 



[G opc {x) - r o °° W(x,q Q )p(q a )dq } 
a 2 (x) 



(25) 



The error a{x) in the above equation is determined from 
the uncertainty of the vacuum condensates and is eval- 
uated by the statistical method explained in (l8j |. S[p] 
is known as the Shannon-Jaynes entropy and can be ob- 
tained as 



S[p] 



dqo[p(qo) ~ m(q ) 



p( go )log(^))], (26) 
m(g ) 



where the function m(oj) is called the default model and is 
an input of the MEM method. It namely corresponds to 
the function which maximizes P[p\GH] if no information 
from the OPE is available (or its uncertainty is infinitely 
large). We therefore use m(oj) to specify the values of 
the spectral function at very high and very low energy, 
for which Gope(2;) provides only a weak constraint. 

Since P[G|i/] does not depend on p(us), it can be re- 
garded as a simple normalization constant and we can 
thus ignore it. From the preceding equations, we now 
have the specific form of P[p\GH]: 

P[p\GH] cx P[G\pH\P[p\H] 

- p QIp] 



where 



Q[p] = aS[p}-L[p]. 



(27) 



Therefore, to get the most probable p(w), we have to 
solve the numerical problem of obtaining the form of p(uj) 
that maximizes Q[p]. 

As for the free parameter a, we treat it in accordance 
with the conventional Bryan algorithm, which means 
that it is integrated out at the end of the calculation 
(see J28| for details). Using this approach, we can get 
the most probable spectral function and additionally its 
error (dp), averaged in some specific energy region. In 
the figures of this paper, we show the errors as three hor- 
izontal lines, whose lengths stand for the region of qo, 



r 


0.5 


1.0 


1.5 


2.0 


Smin of positive parity 


-1.91 


-4.32 


-6.75 


-9.19 


s m in of negative parity 


-1.22 


-3.10 


-5.02 


-6.94 



TABLE II: Values of s min [GeV 2 ] at a 
and fixed r fGeV 4 l. 



1.0 and p 



-0.9 



for which the error (5p) is calculated, while their heights 
correspond to (p) + (Sp), (p), (p) — (Sp), respectively, (p) 
representing the value of the spectral function averaged 
over the region specified by the horizontal length of the 
error bars. For more details about MEM, we refer the 
reader to [13, [H. 



IV. ANALYSIS OF MOCK AND OPE DATA 

Firstly, let us explain how the default model m(qo), 
the analyzed parameter regions of r and s and the values 
of a and /3 are determined. 

As mentioned in the last section, m(qo) should reflect 
our prior knowledge of the spectral function such as the 
asymptotic behavior in the high or low energy region. 
Hence one may naturally conceive three types of the de- 
fault models. The first one is a constant consistent only 
with the asymptotic behavior of the spectral function at 
low energy (which is close to zero as the spectral func- 
tion does not have any strength below the nucleon ground 
state), the second one is also a constant which reflects the 
asymptotic behavior at high energy (which is obtained 
perturbatively) and the third one is a combination of 
the first and second with the correct behavior at both 
high and low energy, which we call as the hybrid default 
model. Carrying out the analyses with three types de- 
fault models, we have confirmed that the position of the 
lowest peak in the extracted spectral function does not 
much depend on this choice. Therefore, since only the 
hybrid model reflects the correct behavior of the spectral 
function at both high and low energy, we will in this pa- 
per only show the results using this default model. Its 
specific form is parametrized as follows: 



TOhybr(w) 



5 + 2/3 + 5/3 2 



1 



128(2tt) 4 1 + e (9th-«o)/5 ' 



(28) 



The concrete values of gth and 8 will be determined in 
the next subsection. 

Next, we discuss the employed parameter region of 
r and s, which is mainly restricted by the convergence 
of the OPE. Since the calculation of the the correlation 
function using the OPE is truncated at a certain order, 
one should keep the contribution of the highest dimen- 
sional term small to justify the truncation and to reduce 
the possibility that higher order terms spoil the result. 
Therefore, we use the well established criterion that the 
ratio of the highest dimensional term is less than 0.1 of 
the whole Gqpe(x). From the property of the kernels, we 
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q th =2.5GeV, 5=0.1GeV q th =3.0GeV, 5=0.1GeV q th =3.0GeV, 5= 0.2GeV 




0123012301234 
q [GeV] q [GeV] q [GeV] 



FIG. 6: Spectral functions extracted from the MEM analyses of G moc k(a;) using various default models at a = 1.0 and /3 = —0.9. 
The solid, dotted and dashed lines show p(qo), the default model and the mock data respectively. The values of parameters 
employed for the default model are given on top of the figure at each row. The parity of the corresponding spectral function is 
shown on the top right corner of each figure. 



expect that the spectral function at small values of r will 
be more sensitive to narrow structures such as the lowest 
peak, while the spectral function at larger values will to 
a large extent be fixed by the continuum. Hence, to ex- 
tract as much information as possible from Gope(^), we 
use several values of r at the same time (r = 0.5, 1, 1.5, 
2 GeV 4 ) and determine the lowest value of s (denoted as 
s min) a t each t. The maximum values of s (denoted as 
Smax), are set to s max = s m i n + 1 GeV 2 . For the detailed 
procedure of the MEM analysis using simultaneously two 
adjustable parameters, we refer the reader to (l9j . 

As for the value of a and j3 of Eq.(JT3J), a can be fixed to 
1 without loss of generality since the interpolating fields 
are determined by the ratio of a and (3. Hence we set 
its value to 1 in the following. With respect to /3, we 
find that the most reliable results can be obtained around 
f3 = —0.9. For other values of (3 such as and oo, the con- 
tribution of the mixed condensate to the OPE becomes 
large and thus the MEM analysis only leads to spectral 
functions consistent with the input default model. This 
happens because of the large uncertainty of the value of 
the condensate, which weakens the constraint of the like- 
lihood function of Eq. (f2"5")) . Near j3 = 1, the subleading 
dimension 6 term proportional to (qq) 2 turns out to be 
the largest contribution to the OPE. Because of the large 



uncertainty related to this term, this is not a very useful 
sum rule. Therefore, we carry out the analyses at a = TO 
and {3 = —0.9. The respective minimum values of s at 
each r determined by the above criterion are given in 
Table. HI 

Finally, let us briefly make a comment on the high 
dimensionality of the spectral function for bayonic chan- 
nels. It is easily understood from dimensional consider- 
ations that unlike in the meson case, the contribution of 
the continuum states to the baryon spectral function is 
proportional to q$ and thus strongly enhanced. As was 
done in similar studies using MEM and lattice QCD, we 
will analyze p(qo)/qQ instead of p(qo). we denote p(<7o)/<7o 
as p(qo) from now on, meaning that Eq. (|22[) changes to: 

/>oo 

Gope{x)= W(x 7 q )q 5 p{q )dq . (29) 
Jo 



A. Analyses using mock data 

In order to clarify the effectiveness of the MEM analy- 
sis and to obtain suitable values for the parameters deter- 
mining the default model, we first analyze the mock data 
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FIG. 7: The positive (left) and negative parity (right) spectral functions extracted from MEM analyses of the OPE data 
Gqpb (s, r) of Eqs. (|18H20l) at a =1 and /3 = —0.9 (solid lines). The dashed lines show the respective results for omitted a s 
corrections. The parity of the corresponding spectral functions is shown on the top right corner of each figure. 



constructed from a phenomenological spectral function of 
a simple "pole + continuum" form: 

Pmock{qo) ~ m1 6{Q °- Ml) + 128(2^ 1 + e ^l- 

(30) 

Here, we use uj = 1.3 GeV and 5 = 0.05GcV. For the 
positive parity (negative) states, the mass of the lowest 
lying state Ml and the residue Xl are taken as Ml = 
0.94 (1.535) GeV, A L = 5.0 x l(r 5 GeV 6 . We apply the 
MEM analysis to G mock (x) = dq W {x , q ) p nlock {q ) 
instead of Gope(x). For the error entering Eq.(|25|l. we 
use the one obtained from the OPE data to mimic the 
actual analysis to be carried out later. The results are 
shown in Fig. [6j In all cases, the peak is successfully re- 
constructed at about the input value Ml, the differences 
between Ml and the position of the obtained peaks being 
less than 40 MeV. On the other hand, the continuum is 
not reproduced well and its behavior strongly depends on 
the values of the parameters of the default model. This 
behavior can be understood by considering the properties 
of W(x,q ). Firstly, one notes that W(x,q ) of Eq. (|T5j) 
oscillates with increasing frequency as qo becomes larger. 
Therefore, in the high qo region, the contribution of the 
continuum is largely reduced due to the fast oscillating 
kernel, which makes it difficult to reproduce its proper- 
ties. The details of the default model however do not 
alter the nature and position of the lowest peak of the 
spectral function, which is our main object of interest. 
Hence, we can conclude that there exist several default 



models which are equally valid for the problem studied 
here. One nevertheless has to make a choice and we will 
for the following analysis of OPE data use q t h = 3.0 GeV 
and 5 = 0.1 GeV for both the positive and negative spec- 
tral functions. 



B. Analyses using OPE data 

Having finished all necessary preparations, we can 
now carry out the analysis using the real OPE data 
Gqpe (s, r), given in Eq. (j2TJ|) . The corresponding spec- 
tral functions are shown in Fig. where the results 
with (without) a s corrections are shown as solid (dashed) 
lines. Comparing the solid and dashed lines, it firstly 
should be noted that the qualitative behavior of the spec- 
tral functions does not change much, which shows that 
the nucleon properties are mainly determined by the non- 
perturbative condensates (which is essentially the chiral 
condensate of dimension 3 here), and not by the pcrtur- 
bative corrections. On the other hand, the lowest peaks 
derived from the full OPE of Eg. (f2T)|) are sharper and 
closer to the experimental values than those including 
only the leading order Wilson coefficients. These findings 
indicate that the inclusion of the a s corrections improves 
the accuracy of the sum rules. Furthermore, we can ob- 
serve that the a s terms cause an increase of the residue 
of the ground state, which is qualitatively consistent with 
the analysis presented in Q. 

Let us now discuss in detail the results of the analy- 
sis using the full Eq. (f2"0")) . As could be expected from the 
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OPE data shown in Fig. [5l the positive and negative par- 
ity spectral functions behave differently and the splitting 
between the two lowest peaks is clearly observed. In the 
positive parity spectral function, peaks are found at 950 
MeV and 1930 MeV. As can be inferred from the error 
bars, the lowest peak which corresponds to the nucleon 
ground state is statistically significant, while the second 
one is not. Furthermore, it is observed that while we 
can extract a clear signal of the ground state, no indica- 
tion for the first excited positive parity state (the Roper, 
N(1440)) is extracted, which most likely means that this 
state only couples weakly to our employed interpolating 
field. 

For negative parity, peaks appear at 1540 MeV and 
1840 MeV. The lowest peak appears close to the ex- 
perimentally observed lowest negative parity state, the 
N(1535). It however has to be emphasized that, because 
of its large width, we can not conclude that the lowest 
peak corresponds only to the N(1535), as it could also 
contain some contribution of the higher lying N(1650) 
state. Related to this issue, let us here make a comment 
on the widths of our obtained peaks. From the MEM 
analysis of mock data, shown in the last section, it can 
be understood that the width of the extracted peaks are 
mainly caused by the limited resolution of the MEM pro- 
cedure and thus should be considered to be an artifact 
of the present method rather than being some physical 
effect. As it was discussed in [2(|, this lack of resolution 
originates from the low sensitivity of the OPE data on 
the peak- width, which can thus constrain the value of the 
width only weakly. Therefore, the peaks corresponding to 
N(1535) and N(1650) likely overlap and could thus both 
contribute to the lowest negative parity peak of Fig. [7] 
We furthermore find a second peak in the negative parity 
spectral function, its existence is however much less cer- 
tain because of the large errors involved. We therefore 
do not speculate about its possible physical interpreta- 
tion here. 



V. SUMMARY AND CONCLUSION 

We have analyzed the nucleon spectral function using 
QCD sum rules and the maximum entropy method. The 
first order a s corrections are taken into account along 
with the parity projection. Analyzing the parity pro- 
jected nucleon sum rule, we have faced several problems. 
Firstly, there were some issues of technical nature in the 
approach proposed in leading to potential ambigu- 
ities in the higher order terms of the OPE. Secondly, 
and most severely, employing conventionally used kernels 
(Borel, Gaussian), we found that the a s corrections tend 
to be too large for the perturbative expansion of the Wil- 
son coefficients to converge, which considerably lowers 
the reliability of the analysis. Moreover, the continuum 
gives a large contribution to these sum rules, which makes 
the extraction of the information on nuclconic states dif- 
ficult. 



To solve these problems, we have firstly improved the 
parity projected sum rules using the forward correlation 
function and have, by this, placed the approach of [l|| 
on a mathematically solid ground. We then have con- 
structed the parity projected sum rule for the nucleon 
including all known first order a s corrections. Further- 
more, by using the phase-rotated Gaussian kernel, we 
could suppress the a s corrections and could thus improve 
the convergence of the corresponding perturbative expan- 
sion and, at the same time, could reduce the continuum 
contribution. Finally, for investigating the spectral func- 
tion, we have used the MEM which does not assume any 
specific functional form for the spectral function, such as 
the "pole + continuum" ansatz. 

In order to check the effectivity of the MEM technique, 
we have first carried out analyses using mock data and 
have confirmed that the lowest peak of the spectral func- 
tion can be reliably extracted without strong dependence 
on the default model (which is, however, not the case for 
the continuum). Next, we have analyzed the OPE data, 
obtained the spectral functions of both positive and neg- 
ative parity and found that both spectral functions con- 
tain significant peaks near the experimental values of the 
lowest lying states. As the term proportional to the di- 
mension 3 quark condensate, which switches its sign de- 
pending on the parity, dominates the OPE used here, this 
sum rule result provides novel evidence for the scenario 
in which the spontaneous breaking of chiral symmetry 
causes the splitting between positive and negative parity 
baryon states. 

Related to the subjects studied in this paper, there 
are, however, still a few unresolved problems left. One 
such issue is the Roper resonance N(1440), which we are 
not able to find in the present analysis. One possible 
reason for our non-observation of the Roper is that in 
this study, we can only use an interpolating field opera- 
tor with /3/a = —0.9 to obtain reliable results. Hence, 
we cannot exclude the possibility that the N(1440) could 
be observed at other values of /3/a. For such an analysis 
to become feasible, it would first of all be necessary to 
determine the dimension 5 mixed condensate and the di- 
mension 6 four-quark condensate with higher precision, 
which could be possible in future lattice QCD studies. 
For the negative parity states, we presently cannot de- 
termine whether the obtained peak corresponds to either 
N(1535) or N(1650) due to the limited resolution of the 
MEM analysis, which is closely related to the error of 
the OPE data which comes from the uncertainty of the 
various vacuum condensates. Therefore, if the values of 
these condensates are determined more accurately in the 
future, a more detailed analysis of the nucleon spectral 
functions of both parities may become possible. 
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Appendix A: Construction of the sum rule from the 
forward correlation function 



In this appendix, we discuss the details of the calcu- 
lations for obtaining general parity projected baryonic 
QCD sum rules. 



1. Derivation of the forward correlation function 

We derive in this subsection the forward correlation 
function from the time ordered correlator. In the course 
of this calculation, we will repeatedly make use of the 
Fourier transform of the Hcaviside step-function: 



^ )= hS dk °k^ 



Ak a x 



(Al) 



Firstly, let us consider the phenomenological side of the 
forward correlator, which can be calculated as follows: 



IL toI (q )=i / <Pxe"> x 0(x )(0\T[T,(x)rj(0)]\0) 



9=0 



c e(x ) 



d 4 P 



dm x 



p+(m)- 



p 



o 

— m 



m^ 



9=0 



1 

p+{m) 



1 



ko - ie jo 
Oo + <7obo - q 



dm 

■ 7 + m 



P-(m) 



(k + q ) 2 - q 2 - m 2 + ie 
(ko + qoho -q-y-m 



(fco + qo) 



q 2 



9=0 



(A2) 



Let us here briefly explain the difference between this 
paper and the previous study of (lfij . We have here re- 
placed ri(x)rj(0) by T[r](x)rj{0)], where T stands for the 
time ordered product. This change is implemented be- 
cause, strictly speaking, it is only the T-product, for 
which perturbative calculations using Feynman diagram 
techniques are allowed. 

The integrand of the above result contains three poles, 
two in the upper half of the imaginary plane of fco and 
one in the lower half. Therefore, closing the contour of 



integration of fco in the lower half and using the Cauchy 
theorem, we pick up the residue of one pole and get 



n tor (<zo) = 

p+{m) 
P- (to) 



1 

2./< 



dmx 
1 



q - yjq 2 + to 2 + ie 
1 



7o 



q ■ 7 



3o 



[70 



yj q 2 + to 2 — ie 
<f ■ 7 + m 



dm 



p+(m) 



1 



p- M 



qo - to + ie 
1 



q Q - m + ie 



\J q 2 + to 2 
(70 + 1) 

(70 - 1) 



9=0 



(A3) 



From this result, it can be understood that the forward 
correlator contains only poles in the positive qo region 
and is analytic for Im qo > 0. Furthermore, by applying 
the projection operator 5(70 ± 1) = P ± to the above 
equation and taking the trace over the spinor indices, we 
can obtain Eq.©: 



1. 



-Im 



Tr^n^o)] 



— Im 

7T 



p±{qo) 



00 

P±(m) — 

o q - to + ie 



dm 



(A4) 



Next, we consider the OPE side of the forward cor- 
relator. For a local interpolating field constructed from 
three light quarks, the time ordered correlation function 
in coordinate space can generally be expressed as 

II(a;) = \C X ^ -j—z ; — + ^ 7—0 

L (x 2 — ie)° [x — ley 



^ x ^ ( x 2 _ i e \2 ^2 _ i c + • • • ] ^ 



x — ie 
i 



(A5) 



' (x 2 - ie) 3 



— ie) 



■CP- 



(n) 

Here, Cx are constants containing condensates with to- 
tal mass dimension n and dimensionless numerical fac- 
tors. In the above equation, the quark masses are taken 
as m u = = 0. This equation is only correct as long as 
we work at leading order in a s for the Wilson coefficients, 
because higher order corrections may involve additional 
logarithmic dependencies on x 2 . The forward correlator 
can be derived from above equation using the Fourier 
transformation and the step function: 



Tl for (q) = i / d 4 xe zqx 0(x o Mx) 



9=0 



(A6) 
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For example, in the dimension 5 term, one encounters 
the following equation: 



d 4 x9(x )e lqx 
dx 6(x y qoi 



(x 2 — ie)" 2 



9=0 



[Xr. X 



ief 



9=0 



(A7) 

First, we calculate the integrals over the spatial angles, 
leading to 

2-7T f 

— / dx 6(x )e iq ° x ° x 

191 J (A8) 



dr 



(r — x + ie) 2 {r + x — ief 



D i\q\r 



<r=o 



where we have used the definition r = \x\. Next comes 
the the integral over r, which can be treated in a standard 
way with the help of the Cauchy theorem. We thus obtain 

dx a 8(x )- 



Jx (qo-\q\) 



Xo — It 

^ [ dk [dxo . 1 . — 

2i J J k -iex a ~ ie 



9=0 
Jx (q +k ) 



(A9) 



At this point, we have dropped \q\, as there is no danger 
that the limit \q] — > leads to a divergence. Making use 
of Eq. (|Al[) for the integral over xo, we reach the final 
result: 

dftQ- = — 7r 2 ln(— go — ie) + 7r 2 ln(oo — ie). 

<- qo ko-ie 

(A10) 

Here, we encounter a divergence in the second term, 
which, however, leads to no relevant contribution to the 
imaginary part of the correlator, which is the only quan- 
tity that is needed for the sum rules. We can therefore 
ignore it. Taking the imaginary part, we obtain the final 
result: 
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d 4 x9{x )e iqx 



(x 2 — ie)" 2 
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The other terms of dimension 0, 3 and 4 can be calculated 
in a similar way. Hence, we here show only the results: 
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(A12) 



(x 2 — ie) 3 



q=0 



We here, as above, have taken the limit |<f| — > 0. Note, 
however that in all the above calculations, this limit can 
be taken only after the integral over r has been carried 
out, as otherwise the factor l/\q\ appearing in Eq. (|A8[) 
can not be properly treated. 

For the other terms of dimension 6, 7 and 8, it is more 
convenient to work in momentum space where the time 
ordered correlator can be given as 



U(q)=[C^q 4 ln(- 
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Here, as in Eq. (|A5|) . Cq contains only condensates and 
numerical factors. Furthermore, we here have neglected 
all polynomials in q 2 , as they are not relevant for the 
present analysis. The forward correlator can now be de- 
rived from the momentum expression using the Fourier 
transformation and the step function as follows: 
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For the dimension 6 term, this form is the same as the 
phenomcnological side of Eq. (IA2l) . Thus, we obtain the 
final result immediately: 
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The other terms containing condensates with dimension 
7 and 8 are transformed into the corresponding terms 
of the Forward correlator through the same procedure. 
The only difference is that due to the larger powers in 
the denominator, the calculations when using the Cauchy 
theorem become more complicated. We here show only 
the corresponding results: 
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As is clear from these expressions, we can at this point 
not take the limit q — > as it would lead to a divergence. 
This problem is only cured after the integral over qo is 
carried out, as will be discussed in the next section. 

Next, we consider the a s correction terms. For the nu- 
cleon sum rules, these corrections have been calculated 
for the Wilson coefficients of dimension 0, 3 and 6. They 
partly have the same analytic structure as Eqs. (|A5P and 
flABJ) , which can be treated as explained above (we ig- 
nore the running of a s here). However, for the dimension 
and 6 terms, the a s corrections involve contributions 
with additional logarithms, which complicate the calcu- 
lation. In momentum space, these terms are expressed 
as 
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For the term of dimension 0, we use the equations given 
below M : 
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where 7 is the Euler constant. Making use of these two 
equations, we first derive the a s correction term in coor- 
dinate space. Then, we can calculate the forward corre- 
lator using Eq. ([A12p and 
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Note that for obtaining Eq. (|A19[) with the help of the 
Cauchy theorem, one must consider the complicated con- 
tour avoiding the branch cut of the logarithmic function. 

As for the dimension 6 term, we can derive its forward 



correlator as follows: 



Im 



d A xe{x )e lqx 



d 4 p 



(2tt) 4 



-ipx 



jr + ie 



ln(— q — ie) 
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70 Im 



dq l_ 

2iri qo — ie 



qo + fco 



(<7o + h ~\q\+ i() {qo + ko + \q\ - ie) 
x (ln(\q\ - q - k - ie) + ln(\q\ + q + k - ie)j 



q=0 

(A20) 



In the first (second) term, there is a branch cut due to the 
logarithmic function below (above) the real axis. There- 
fore, we can choose to close the contour for each term 
such that it does not contain the cut and then simply 
use the Cauchy theorem. We thus obtain 



Im 



d 4 x9{x )e lqx 
x ln(-q 2 - ie) 



d 4 P _-i 



(27T) 4 

= 7o Im 



p 2 + ie 



1 



9=0 



qo + if 



ln(-<7o - ie) 



(A21) 



In summary, we show the relation between the time or- 
dered correlator and the imaginary part of the forward 
correlator in the first two columns of Table IIIII 

From the above procedure, we can obtain the forward 
correlator of the nucleon channel up to dimension 8 and 
first order a s corrections from the time ordered correlator 
0,[H as: 



U iol (q ) = 7o nr (qo) + nr (qo) 



for , 



jfor , 



Im 



n[° r (qo) 



g=0 
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Im 
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= (Co + C 0a3 )q 5 6(q ) 
+ C 4 ( — G 2 )qo 9(qo) 

IT 

+ C 6 (qq) 2 5(q ) + C 6as (qq) 2 
+ C 8 (qq)(qgcr ■ Gq), 

(C 3 + C 3a J(qq)q 2 9(q ) 

+ C 5 (qga-Gq) 8(q ) 



Cr(qq)(^G 2 ) 

7T 



(A22) 
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with 
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Cs = -* 2^ X 
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Im 



1 



1 
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\q\-ie (q - \q\ + ie) 2 
where /i is a renormalization constant. 

2. Construction of the sum rule 



(A23) 



results obtained above. Here, we assume a general form 
of the kernel, which we take to be either an even (odd) 
function of q and denote it as W e (qo) (<7oWe(<Zo))> wnere 
W e (g , g) is analytic in the upper half of the imaginary 
plane and real on the real axis. Furthermore, we will 
prove that the potential divergences of the higher order 
terms of dimension 7 and 8 vanish through the integral 
of Eq.(Hnj). 

Firstly, we consider the terms up to dimension 5. For 
these, the imaginary parts of the forward correlator is 
proportional to a step function and one therefore just 
has to substitute our findings of the last subsection into 
Eq. ([TU)) and evaluate the integral from to oo. This leads 
to the results presented in Eq. (f20"|) . 

Next, let us discuss the more complicated terms of di- 
mension 6, 7 and 8. As a first step, we show that, if one 
takes the limit \q\ 0, the imaginary parts of the lead- 
ing terms of dimensions 6,8,... (7,9,...) are even (odd) 
functions of qo- By following the procedure of the last 
subsection, it is noticed that all terms appearing at di- 
mensions 6,8,... can generally be written down as 



1 



1 



k + qo 



ti J °'k — ie [(k + qo) 2 + ie]' 
1 f .. 1 ko + qo 



Fi(q ) 



dko- 



2iri J " k - ie (ko + qo + ie) n (k + qo - ie) n ' 

(n=l,2,...) 

(A24) 

Here, we are ignoring any proportional real constant, in- 
cluding 7q. Similarly, for dimensions 7,9,..., we get 



7T- / dk - 



2ni J 'k - ie [(k + q ) 2 + ie}' 1 



F 2 (qo) 
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1 



2-Ki J ko — ie (k a + qo + ie) n (k + qo - ie) n 

(n = l,2,...) 

(A25) 



In this subsection, we show how the final form of the 
sum rule is constructed, starting from Eq. (|10[) and the 



Next, we take the imaginary parts and, after some simple 
manipulations, get for F\(qo) 



ImFi(go) = - FM - F x (q Q ) 
i 



1 



dkr 



1 



1 



2ir J " V ko + ie ko — ie 
while the result for ^2(90) is 
1 



k + qo 
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(A26) 



(ko + q + ie) n (ko + q Q - ie) n (-k a + qo + ie) n (-k + q - ie) n - 



lmF 2 (g ) = - F 2 (qo) ~ F 2 (q ) 

= -7T I dk ^(TrT- + T~~~ 

2ir J V ko + ie ko — ie 



1 



1 



(A27) 



(k + qo + ie) n (k + q - ie) n {-ko + qo + ie) n (-k + q - ie) r 
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Having the above equations at hand, it is now a trivial 
matter to show that 



and 



ImFi(-qo) = -IraFi(go)) 



lmF 2 (-q ) = ImF 2 (g ) 



(A28) 



(A29) 



which proofs that the imaginary parts of the terms of 
dimension 6,8,... (7,9,...) in the OPE of the forward 
correlator are even (odd) functions of go m the limit \q\ — > 
0. 

Therefore, we can immediately evaluate half of the con- 
tributions of the these terms in Eq.(|10j). Namely, we find 
that the leading terms of dimension 7,9,... vanish if the 
kernel is an even function while the terms of dimension 
6,8,10,... vanish if it is an odd function. 

Among the remaining terms, the leading dimension 6 
term with an even kernel can be obtained trivially, as the 
the imaginary part of the forward correlator is propor- 
tional to a delta function at the origin. To calculate the 
non-vanishing part of dimension 7, we consider the case 
where W(qo) is an odd function: 



Im 



d 1o 7T 



1 



2\q\-ieq - \q\ + it 



q W e (q 2 ) 



= - g We(O). 



(A30) 



It can be understood here that the seemingly divergent 
part 1/q disappears after the final integration over go. 
For dimension 8, we consider the case where We(go) is 
an even function: 



Im 



, 1 1 
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Im 
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,11 1 d 
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We(q 2 ) 
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7T d 



2 d(q 2 ) 



We(0). 



(A31) 

After integrating by parts, which leads to the second line 



in the above equation, we see that, as for the dimension 
7 case, the diverging term vanishes. We have explicitly 
checked that for the dimension 9 term, all potential di- 
vergences vanish in a similar fashion. Thus, the same 
procedure could presumably be continued to even higher 
orders, but this is not of much practical use as the higher 
order terms have large numerical uncertainties from the 
corresponding condensates at present. 

Finally, let us discuss the part of the a s correction 
term of dimension 6, which contains a logarithm. Due 
to this logarithm, the contributions of this term using 
both W e {qo) and goIE e (go) are finite since this term is 
neither an even nor an odd function of go- For W e {(^), 
the contribution can be calculated with the help of the 
Cauchy theorem by considering the contour avoiding the 
pole. The result can be obtained as follows: 



Im 



oo 

dqo — 

-oo Qo + it 



ln(-go-ze)W e (g 2 ) 



= -7rln(e)W(0-7r / dq -W e (q 2 ) (A32) 

Je QO 

. [°° , . , , dW e (q 2 ) 
= 2ir dg g lng ) — — 2 — . 
Jo dqfi 

Note that in the second line, t stands for the infinitesi- 
mal radius of the path around the origin of the imaginary 
plane and that the diverging parts of the two terms ap- 
pearing there exactly cancel, leaving only the finite result 
of the third line. 

For qgWe^q 2 .), the contribution can be immediately ob- 
tained as: 



Im 



dqo 



1 



go + it 



In (-g - it)q W e (q 2 ) 



dqo 0(q )W e (q 2 ) 
dq W e (q 2 ). 



(A33) 



-oo 
oo 



All the results of this subsection, using both even and odd 
kernels, are summarized in the third and fourth column 
of Table Hill 
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